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Abstract 

The dynamics of a baby skyrmion configuration, in a 
model Landau-Lifshitz equation, was studied in the presence of var- 
ious potential obstructions. The baby skyrmion configuration was 
constructed from two Q = 1 hedgehog solutions to the baby skyrme 
model in (2+1) dimensions. The potential obstructions were created 
by introducing a new term into the Lagrangian which resulted in a 
localised inhomogeneity in the potential terms' coefficient. In the bar- 
rier system, the normal circular path was deformed as the skyrmions 
traversed the barrier. During the same period, it was seen that the 
skyrmions sped up as they went over the barrier. For critical values 
of the barrier height and width, the skyrmions were no longer bound 
and were free to separate. In the case of a potential hole, the baby 
skyrmions no longer formed a bound state and moved asymptotically 
along the axis of the hole. It is shown how to modify the definition of 
the angular momentum to include the effects of the obstructions, so 
that it is conserved. 
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I Introduction 



The scattering of particles off potential holes and barriers in classical and 
quantum mechanical systems are seemingly different. In a classical system, 
if the particle has sufficient energy it can traverse a barrier; if it does not, it 
gets reflected. When it encounters a hole it speeds up as it passes over the hole 
and is always transmitted. In both cases the particle is either transmitted 
or reflected. In quantum mechanical scattering the particle can be reflected 
and transmitted for either a barrier or a hole but these events occur with a 
certain probability which is dependent on the particle's energy and on the 
size of the barrier or hole. 

In this paper we examine the scattering properties of a topological soliton 
in magnetic systems, whose motion is then governed by the Landau-Lifshitz 
equation. Topological solitons are, of course, classical objects (as they sat- 
isfy classical equations of motions). However, they describe extended objects 
and, as shown in [TJ, some of their properties resemble those of quantum 
systems. Hence in this paper we study this problem further; this time con- 
centrating on systems whose dynamics is described by a Landau Lifshitz 
equation. This equation arises in the dynamics of magnetic bubbles and so, 
in this paper, we look at the behaviour of topological solitons in the presence 
of potential obstructions. The topological solitons under investigation are 
baby skyrmions which are thought to describe the experimentally observable 
magnetic bubbles. 

Our investigation should also shed some light on the properties of mag- 
netic bubbles. Such bubbles are not mathematical artifacts but have been 
produced experimentally by subjecting a ferromagnetic material to a pulsed 
magnetic field. The ferromagnetic domains of the material are then squeezed 
by the field. If a field of large enough magnitude is applied over a sufficient 
length of time the magnetic domains of the system collapse leaving behind 
the material which is uniformly magnetised in the direction of the applied 
field. This process is not instantaneous or uniform. As the domains tend 
to align with the applied field, the pulsing of certain materials results in the 
pinching of the domain walls into a cylindrical domain called a magnetic 
bubble. Magnetic bubbles were once considered as an alternative form of 
memory storage due to the density at which they could be stored. 

In real three dimensional systems the magnetic bubbles are stabilised by 
the finite thickness of the thin films in which they are created. The model 
examined here is strictly in 2-D. In this case the bubbles are stabilised by 
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the introduction of the slightly artificial Skyrme term. For more realistic 
calculations one would need to perform 3-D simulations. Magnetic bubbles 
and their properties have been extensively experimentally researched; further 
information on this research can be found in [3] and [3]. 

The baby skyrme model in (2+l)-dimensions is defined by: 



where is a 3-component scalar field and the indices run over the space-time 
coordinates. This is referred to as the baby skyrme model to distinguish 
it from the full skyrme nuclear model of baryons. The first term is the 
exchange energy, the second and third terms are the skyrme term and the 
potential term respectively. The latter two terms have been introduced to 
avoid the consequences of Derrick's theorem [5] and to stabilise topological 
soliton solutions in 2-dimensions. The condition that <fi 2 = 1 is imposed so 
that the target space is the 2-sphere, such that <fi is now a map : M 2 — > S 2 . 
For finite energy solutions it is necessary for the fields to tend to a vacuum 
at infinity, where 03 = 1 at oo. This results in a compactification of IR 2 so 
that (J) now takes values in the extended plane M 2 U oo, which is topologically 
equivalent to S* 2 . The constraint equation 2 = 1 and the boundary condition 
at infinity results in the field <fi becoming a non-trivial map : S 2 — > S 2 . 
Each soliton solution is grouped into a different homotopy class according 
to the winding number, or topological charge, of this map. The topological 
charge Q is given by: 



where the indices i, j run over the space coordinates and Q G Z. The topolog- 
ical soliton solutions of the baby skyrme model are called baby skyrmions. 
Here, for simplicity, we shall refer to baby skyrmions of charge Q as Q- 
skyrmions. Recent work has shown that the continuum dynamical equation 
in anti-ferromagnetic systems, also resembles a second order relativistic wave 
equation [2] in which such solitons solutions exist. The work in [2] established 
many interesting properties of solitons in such systems. 

Our discussion so far has concerned the relativistic systems of skyrmions. 
However, as we said above, they also arise in the description of magnetic 
bubbles, but this time their evolution is described by the first order Landau- 
Lifshitz equation. The Landau-Lifshitz equation is given by: 
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where W is the energy functional written as: 

p POO 

W — wdxdy, (4) 



and w is the static part of ([I]) given by: 

w = ^0 • <9,0 + \li\{di<t ■ d 4f ~ ( d 4 ■ d M94 ■ dM + V(4>). 



Thus we can write down 



sw . 



sw a dv(4>) 

+ i 72 {29 4 [(9,0 ■ ^0)^0] - $ [(^0 ■ 9,0)^0] - ^ [(dtf ■ 9,0)^0] }. 

Analysis of the dynamics in Landau-Lifshitz systems has been greatly 
simplified by the work of Papanicolaou and Tomaras [TU], who constructed 
unambiguous conservation laws for the system governed by (j3J) • In their work 
they found that the important quantity was the topological charge density 
q: 

q = e^-0 • (<9j0 x <3i0). (5) 
The conservation laws are constructed from the moments of q. They involve: 



oo 



I = - 1 1 x 2 qdxdy, (6) 



oo 



oo 



m = \\ (<f) 3 -l)dxdy, (7) 

J J —oo 

J = l + m, (8) 

where / is the orbital angular momentum, m is the total magnetization in 
the third direction and J is the total angular momentum. Conservation laws 
for the system were constructed by examining the time evolution of q: 

q = - eij d.d l(T . h (9) 
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where d\Oji can be written in terms of the energy functional W: 

fsw \ 

d ^oi= TX-W' ( 10 ) 



Taking an explicit time derivative of ([6]) gives: 



2 

and we note that this can be recast as 



/' = - / / x 2 qdxdy, (11) 



1 = > ,jV,j dxdy, (12) 



oo 



by integrating (fTTl) by parts. In the case of a system with a symmetric tensor 
(Jij ,1 = and angular momentum is conserved. 

The guiding centre coordinate R of the soliton is defined as the first 
moment of the topological charge density q: 



R=- — — // xqdxdy. (13) 

Since our solitons are in two spatial dimensions, their position can generally 
be defined as the location of the centre of each soliton. This definition can be 
interpreted in two ways. One can consider the soliton centre to be the point 
at which the third component of the field 03 = —1 or to be the maxima of 
the topological charge density q. It was seen in previous simulations [8] that 
both of these definitions produce near identical trajectories. One can also 
consider the mean squared radius r of the solitons, defined by: 



1 P /-oo 



One can expand out (Tl4"|) to find a relationship between the mean squared 
radius of the solitons and /: 



/• 



2 * D 2 



2ixQ 



-R 2 . (15) 



This relationship between I, r and R greatly helps to understands the 
dynamics described in later sections. 
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Much of the previous work on baby skyrme models, in the context of Landau- 
Lifshitz dynamics, has concerned the choice of two different potential terms 



Models which employ either of these potentials are commonly referred to 
as the holomorphic baby skyrme model (fl6j) and the 'new' baby skyrme model 
(I17jl to distinguish it from the holomorphic one. The holomorphic model 
was first studied in the context of Landau-Lifshitz dynamics in [8], since it 
provided an analytical solution to the system of equations. The topological 
solitons of this model are polynomially localised. In [8] it was shown that two 
1-skyrmions orbited around each other along deformed circular trajectories. 
Their work involved a local magneto-static field in addition to the three 
previous terms in (TjQ), which resulted in the non-conservation of / and m. 
The total angular momentum J was well conserved in time. The authors of 
[8] attributed the non-conservation of / and m to the non-symmetric structure 
of dji due to the presence of the magneto-static field. 

The new baby skyrme model is a more realistic case of easy axis anisotropy 
and we use this in our study i.e. we consider V((f>) to take the form of (fiT|) . 
Some of the work done on it has been in relation to the dynamics of magnetic 
bubbles [7J; this study also involved a local magneto-static field. Analytic 
solutions do not exist and so all solutions must be found numerically. The 
topological solitons of this model are exponentially localised. In [7J it was 
found that two Q = 1-skyrmions orbited each other on a circular trajectory 
modified by a Larmor procession due to the magnetic field. J was well 
conserved in time but its constituent components / and m were not. The 
arguments for the non-conservation were identical to those in the holomorphic 
model. Since those early papers, there has been extensive work done on both 
these models and their multiskyrmion structures. The details can be found 
in the work of Weidig [6], or for a larger class of potentials V(<f>) in [9]. 

II Constructing the initial field configuration 

Seeking a static field configuration which is a solution of ([3]) for a poten- 
tial term V(4>) of the form ffTTl) . we assume the solitons take the form of a 



V {</>): 




(16) 



(17) 
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hedgehog configuration for a Q-skyrmion: 



= (cos(Qfl) sin(/(r)), sin(Q#) sm(/(r)), cos(/(r))), (18) 

where f(r) is the profile function satisfying certain boundary conditions, 
9 is the polar angle and Q is the topological charge of the skyrmions. The 
skyrmion solutions are minima of the energy functional (J4]) . We are interested 
in the dynamics of Q = 1 skyrmions. Inserting the hedgehog configuration 
for Q — 1 into the energy functional and minimising the integral, results in 
a second order differential equation for the profile function f(r): 

( 72 sin 2 A / 72 sin 2 / 
/ I 7ir + : ) + / I 7i ~2 — 



, 2 /'72sin/cosA 7!sin/cos/ 
+ / 7 3 r cos / sm / = 0. (19) 



This can be rearranged into the form /" = h(r,f,f). The profile function 
must satisfy certain boundary conditions for there to exist finite energy so- 
lutions. The boundary conditions impose constraints on f(r) at the origin 
and infinity: /(0) = 7r and /(oo) = 0. The second order differential equation 
for fir) can be solved numerically using the shooting method. With the 
profile function obtained we can construct from (fT8j) a 1-skyrmion solution 
to (J3J). A two 1-skyrmion configuration can be constructed by the superpo- 
sition procedure. The easiest method to do this is to transform the fields to 
a stereographic variable Q: 



n 



(pi + i<t>2 

1 + 03 ' 

«i + |n| 2 
i- \n\ 2 



i + \n\ 2 

The stereographic variable Q can then be rewritten in terms of the hedgehog 
anzatz variables, f(r) and 9 as: 

Q = tan(^)e ie . 
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To construct a two Q = 1-skyrmion configuration in which the two skyrmions 
are in an attractive channel, we take: 



Q = fix-fta, (21) 
ni = tan(^)e iei , 



n 2 = tan (lM) e ie^ 



where r, = yj (Jx — Xi) 2 + (y — Hi) 2 ) and 6i = Tan 1 (fzf 1 ) are calculated 

relative to the centres of the skyrmions (xj, yi). During the simulations it was 
found that the configuration constructed in this manner did not replicate the 
skyrmion ring configurations for small values of 



d=J ((xi — X2) 2 + (yi — y2) 2 )] see [6]. The superposition procedure was a 
very good approximation to a two 1-skyrmion configuration for values of d > 
6, where the skyrmions were well separated to be distinct. To obtain a true 
representation of the configuration for all values of d, we used a gradient flow 
method to 'relax' the field configuration. The above field configuration (|2~TT) 
constructed for skyrmion separation d = 8, was used as an initial condition 
of the gradient flow equation given by: 

d(j) SW 7 t 

Tt = ~ K ^l + k ^ 

where k is a Lagrange multiplier introduced such that the constraint 2 = 1 
is satisfied. The field configurations obtained by this relaxation method show 
the required ring like properties for small values of the skyrmion separation. 



Ill Potential obstruction 

In this paper we study the scattering properties of a two Q = 1-skyrmion 
configuration on a potential obstruction which is localised in a finite region 
of space. In constructing the obstruction we adopt a similar approach used 
in the previous work of one of the authors [1] and introduce a term into the 
Lagrangian (CQ) which vanishes in the vacuum state 03 = +1. The obstruc- 
tions need to be introduced in this way so that the tails of the solitons are 
not changed by the obstruction. Therefore, we add the additional potential 
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term Vobstructioni^z) which is identical to the potential in ([T]) and so it effects 
the Lagrangian of ([I]) by changing the potential coefficient from a constant 
to one that depends on the space coordinates. The new Lagrangian of the 
system written in terms of the previous one: 

l~-new l~-old Vobstruetion\4 l 3) 5 

where: 

V bstruction((p) = T^i 1 ~ 4>l)- ( 22 ) 

T can be either positive or negative. The sign of V determines whether the 
potential obstruction is a hole or a barrier. When T > 0, the obstruction is a 
barrier. Conversely when T < the obstruction is a hole. The introduction 
of this term, implies that (CD) can be rewritten as: 

£ = -^Ixd4-d4--^ 2 [{d4 ■ dif) 2 - (dit-djfyidit-djcjj)] - TpO, 

(23) 

where the potential term coefficient 73 is now a function of the coordinates 
(x, y) and the static part of (pQ) is only considered as imposed by the Landau- 
Lifshitz equation ([3]). This inhomogeneity will be localised to a finite region 
of space. The value of 73(2:,?/), in this region, will determine whether the 
obstruction is a hole or a barrier. In this region, 73(0:, y), can therefore be 
summarised as: 



T3V*^> V^region 73 T, 



T > Barrier, 
T < Hole. 



IV Numerical procedures and the free sys- 
tem dynamics 

Unfortunately, it is impossible to solve ([3]) analytically we have therefore 
had to study this problem numerically. The fields and their derivatives were 
discretised in the usual manner and were placed on a lattice of 251 x 251 
points, with lattice spacing dx = 0.1. The numerical integration of the 
3-coupled differential equations of ([3D involved the use of a 4 th order Runge- 
Kutta method of simulating time evolution with a time step of dt = 0.001. 
The various integrals calculated throughout the simulations were performed 
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using a 2-D Simpson's rule. The constraint equation requires that the fields 
lie on the 2-sphere, =1, and this was imposed at every time step by 
rescaling each field component so that d>i — > p= . 

The skyrmions were initially placed at (0, ±d/2) in the upper and lower 
planes, where d is the distance between the two skyrmion centres (xi,yi). 
The trajectory of each skyrmion was tracked by following the maxima of 
the topological charge density and interpolating between the lattice points. 
All of the simulations have been performed for a skyrmion separation of 
d = 6. This was found to be the optimum distance, where the skyrmions are 
separated enough from each other to be distinct but close enough to interact. 
The coefficients 7; have been set to unity in all the simulations unless stated 
otherwise. 



IV. 1 No obstruction dynamics 

Initially we examined the behaviour of the skyrmions without an obstruction 
i.e. with T = 0. Figure ([1]) shows the trajectory of the upper skyrmion 
of a two 1-Skyrmion configuration for T = 0. The skyrmions orbit around 
the configuration's centre (0,0). The trajectories of each skyrmion lie along 
a circle of radius r ~ 3. Their position undergoes mild oscillations during 
the simulation. The total energy E tot and angular momentum J = I + m 
are conserved with time. Additionally, each individual angular momentum 
component I and m is also conserved. The time scale for one period is 700 
sees. This motion of two baby-skyrmions in a Landau-Lifshitz system is well 
understood and an analogy with the Hall motion of two interacting electrons 
is usually invoked when discussing their trajectory [TO] . 



V Simple obstruction 

There are many choices one can make for the geometry of the potential ob- 
struction. The simplest choice initially studied was a symmetric obstruction. 
The obstruction was placed symmetrically along the x-axis with width b, i.e. 
starting at y = —b/2 and continuing up to y = b/2 and extending for all 
values of x. In our study we look at the differences of dynamics due to holes 
and barriers. They will be discussed in the following sections. 
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Figure 1: Trajectory of the upper skyrmion in the absence of any obstruction 
for a two 1-Skyrmion configuration. 

V.l Potential hole 

We start by recalling that in the absence of all obstructions the skyrmion 
execute a circular path around their centre. Fig. @ shows the trajectories 
of the upper and lower skyrmion encountering a potential hole for different 
values of T, for b = 1. In all plots the skyrmions initially try to execute the 
trajectory of Fig. ([I]) but are deflected. They move asymptotically along the 
axis of the hole at an approximately constant value, y m ax- It is clear from 
Fig.(T5]) that the larger the |T|, the larger the value of y ma x- The skyrmions of 
Fig.([2])a are able to get 'closer' to the hole than the skyrmions of Fig.(T2])b or 
Fig.(T2])c. It is reasonable to suggest that the repulsion of the skyrmions by 
the hole increases with increasing Fig. ([3]) shows the trajectories of only 
the upper skyrmions interacting with the potential hole for various values of 
r when b=3. Comparing this plot with those in Fig. (|2j) one sees the effect of 
a larger b. The larger the b the larger the y max for a given V. The skyrmions' 
tail can feel the hole earlier in a system of larger b, than in a system of smaller 
b. The same dependence of y ma x on T is evident in Fig. (j3J). In Fig. (j3J) the 
skyrmion trajectories of T = —0.25, —0.5 are seen to be reflected by the 
boundary. The skyrmions in the lower plane execute similar trajectories. 

To explain the observed behavior of skyrmions when encountering a hole, 
which may at first sight, appear as non-classical, we need to examine the 
binding energies of the configuration. Table (jTJ presents the energies of a 
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single skyrmion in ([3]), with 73 set at the same value throughout the whole 
lattice. If the energy of the two 1-skyrmion configuration in the presence 
of a hole is denoted by E 2 and E\ is the energy of a single 1-skyrmion in a 
system with 73 set at the same value as the hole in E2, then the binding, or 
the interaction energy, is given by: 

E B = E 2 - 2El (24) 

In the system with b = 1, T = —0.1 we have E2 = 2.1206/87T and E\ = 
1. 0454/87T. The binding energy of the two skyrmions in this system is thus 
Eb = 0.0356/87T, which is positive. Hence in the case corresponding to our 
simulations the skyrmions are no longer bound. Correspondingly, the system 
with b = 1, T = -0.5 has E 2 = 2.1205/8tt and E l = 0.9494/8vr. The binding 
energy of the two skyrmions in this system is Eb = 0.2217/87T. This is also 
positive but larger by a factor ~ 6.3. Thus both skyrmion configurations 
placed in the presence of a potential hole are no longer bound. This analysis 
also explains some of the other features of the system. The trajectories 
of Fig.Qa take substantially more time than Fig.Qb or Fig.Qc. Eb is 
very much larger when T = —0.5 than for T = —0.1 and, hence, this extra 
available energy allows the skyrmions in V = —0.5 to separate more quickly. 
In consequence, the skyrmions for the systems of larger |T| have a larger 
speed in this 'asymptotic' state than those of smaller 
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Figure 2: Trajectories for both the upper and lower skyrmions with a poten- 
tial hole, for b = 1 and with the time length shown in brackets: a)T = —0.1 
(690 sees), b)r = -0.25 (290 sees), c)r = -0.5 (215 sees). 
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Figure 3: Plot of the trajectories of the upper skyrmion in a system with a 
potential hole of width b=3 for various values of T. 
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Table 1: Table showing the variation in total energy, in units of 8ir, for a 
1-skyrmion with the 73 coefficient set at the corresponding value all over the 
whole lattice. 
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V.2 Potential barrier 



Next we have studied several cases of the scattering of the same two 1- 
skyrmion configuration off potential barriers. Fig.(j3]) shows the trajectories 
of the upper and lower skyrmion of a two 1-skyrmion configuration scattering 
off a potential barrier of width b = 2, for various values of T. It can be seen in 
each plot that the skyrmions are deflected as they traverse the barrier. This 
deflection always occurs in the direction of the centre of the configuration 
and hence the normal circular path is deformed as the skyrmions overcome 
the barrier. Trajectories for a smaller barrier width show a sharper deflection 
than those with larger b. Since the skyrmions are extended objects, when 
they are traversing the barrier they feel the barrier the most at the middle 
point. It is then only natural that the maximum point of deviation from 
the normal circular path will be at this point. Comparing trajectories of 
the same T but with different values of b, we note that the skyrmions in the 
system with the larger value of b have more time to adjust to the barrier 
once they are 'on' top of it and therefore their path is not as sharp as for 
a smaller b. The larger value of b 'smooths' out the sharpening effects seen 
in the system with a smaller value of 6's. More interestingly, during this 
deviation the skyrmions speed up as they traverse the barrier. In Figs.(j4])a 
and Figs.(j4])b the times taken for the skyrmion centre to reach the edge of 
the barrier are 110 sees and 140 sees respectively, but the times taken for the 
centre to traverse the full width of the barrier are only 25 sees and 12.5 sees. 

The deflection and speeding up of the skyrmions can be explained by ex- 
amining Table(T5]) which shows a table of the energies of the two 1-skyrmion 
configurations for differing values of the skyrmion separation d, with the 
potential coefficient 73 set the corresponding value over all the lattice. As 
explained in section(2), our two 1-skyrmion configuration was constructed in 
such a way that the skyrmions were in an attractive channel. When they 
encounter a barrier the energy of their configuration would have to increase. 
The skyrmions counteract this increase due to the barrier by reducing their 
separation distance d. This is clear by considering the energies of the config- 
uration away from the barrier with d = 6. The energy of such a configuration 
is E = 2.1277. If the same configuration was then placed on the barrier with 
T = 0.2, the energy becomes E = 2.1283. Thus any increase in the potential 
energy due to the barrier, must be compensated by a reduction in d. The 
energy increase of the system would be at its greatest when the skyrmions 
are in the middle of the barrier hence the biggest deflection is seen at this 
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point. Due to this quick adjustment of d, the skyrmions speed up as they 
traverse the barrier. 

Another interesting feature of the scattering on a barrier is the 'transition 
dynamics' shown in Fig (jSJ). These plots show the transition to a state in 
which the skyrmions do not traverse the barrier and, instead, move away from 
each other. The transition to such a state is shown through the variation 
in the potential coefficient T for a fixed value of the barrier width 6 = 3. 
Similar plots could also have been obtained by choosing a fixed value of 
T ~ 0.25 and increasing the barrier width b from b = 2 to b = 3. This 
effect is due to the binding energies of the skyrmion configuration. Using 
the previous definition of the binding energy (1241) and its constituent parts, 
one can examine the binding energies in the barrier system. In a system 
with b = 2 and T = 0.1, E 2 = 2.1317/8tt and E x = 1.0697/8*- therefore 
Eb = — 0.0077/87T so the skyrmions are still bound. In b = 2 and T = 0.3, 
E 2 = 2.1397/8vr and E 1 = 1.0702/8vr therefore E B = -0.0007/8tt and the 
skyrmions are still bound although a bit more loosely than for the smaller 
value of T. Next, consider the state where the skyrmions separate from each 
other i.e. for b = 3 and T = 0.25. Then E 2 = 2.1482/8vr and E x = 1.0701/8vr 
and therefore, Eb = 0.0079/87T. Thus, in this system the skyrmions are no 
longer bound and separate from each other. 
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a 




Figure 4: Trajectories of upper and lower skyrmions for b=2 barrier system 
for various values of T: a)T = 0.1, b)T = 0.25 , c)T = 0.3. 
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a 



b 




Figure 5: Plots showing the transition to asymptotic state b=3 T = 
0.25,through variations in T for a fixed b=3: a)T = 0.10, b)T = 0.15, 
c)r = 0.20, d)r = 0.25 . 
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Table 2: Table showing the variation in total energy, in units of 87r, for a 
two 1-skyrmion configuration with the 73 coefficient set at the corresponding 
value all over the lattice. 



VI Angular momentum 

In this section we present the explanation of our results based on the study 
of the total angular momentum J. The orbital angular momentum I and the 
total magnetization in the third direction m were calculated through all the 
simulations using the definitions given by ([6]) and (J7|). Fig.(j6]) shows a plot 
of I, m and J for a two 1-skyrmion configuration interacting with a potential 
barrier of width b = 2 and T = 0.25. It is clear that rh = throughout, but 
that I, and J, are not conserved in time. 

In systems involving a two 1-skyrmion configuration the guiding centre 
coordinate R , defined in (1141) . corresponds to the centre of the configuration. 
A calculation of R during the barrier simulations has indeed shown that 
R = 0. This is expected since the trajectories of the skyrmions in the system 
are always symmetric with respect to a reflection through the origin and so 
the centre of the configuration always lies at this point. Considering (fTST) 
with R = 0, we note that the orbital angular momentum I and the average 
size of the skyrmions r are now directly related to each other in the barrier 
system by: 
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Figure 6: Plots of the orbital angular momentum I, total magnetization m 
and total angular momentum J = I + m for b=2 T = 0.25. 

FigflZj) shows a plot of the average skyrmion radius as a function of time for 
a potential barrier system of width b = 2 and T = 0.25. The points at which 
r(t) approaches its minimum, corresponds to the skyrmions traversing the 
barrier. The point at which they have reached the maximum of the barrier 
corresponds to the minimum of r(t). Thus as the skyrmions traverse the 
barrier their average size decreases from its starting value by around %20. 
Since the tail of the skyrmion is exponentially localised and this localisation 
is governed by the potential coefficient parameter 73, it is expected that due 
to the inhomogenaity in 73 their size would decrease in the region of larger 
73 explaining the observed behaviour in I. 

Fig(JED shows a plot of I, m and J for a two 1-skyrmion configuration 
interacting with a potential hole of width b = 2 and T = 0.25. Again, it 
is clear that rh = throughout, but /, and J, are not conserved in time, 
analogous to what was seen in the system involving a potential barrier. The 
guiding centre coordinate R for this system can also be shown to vanish and 
thus (125j) is valid also in systems with potential holes. Using this, we can 
therefore plot r(t) for a potential hole. FigQ shows a plot of the average 
skyrmion radius as a function of time for a potential hole system of width 
6 = 2 and T = —0.25. It is clear from the plot that as the skyrmions approach 
the boundary asymptotically along the edge of the hole, the average size of 
the skyrmion increases continually, increasing to approximately 3 times its 
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Figure 7: Plot of the average skyrmion radius as a function of time for a 
potential barrier with b = 2 and T = 0.25. 

initial size. This is due to the tail of the skyrmions. The skyrmions cannot 
penetrate the hole, as explained in the previous sections, but its tail can. The 
exponential localisation of the skyrmions, as explained earlier, is governed by 
the potential coefficient 73. In the region of reduced 73, their average size is 
able to grow and so it continues to increase until they reach the boundary 
of the system, where, they get reflected. In our simulations we saw that 
following this reflection the skyrmions' size decreases back to its starting 
value as the system tends to its starting point. 
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Figure 8: Plots of the orbital angular momentum I, total magnetization 
m and total angular momentum J = I + m for a potential hole with b=2 
r = -0.25. 



Figure 9: Plot of the average skyrmion radius as a function of time for a 
potential hole with 6 = 2 and V = —0.25. 
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VII Further discussion of Angular Momen- 
tum 



Here we discuss further the apparent non-conservation of J. Assuming the 
definitions of I, m and hence of J to be valid and shown to be true in the 
free system, we ask ourselves can we explain this more qualitatively? Let 
us consider the behaviour of J. The form of l(t) from Fig. ([6]) indicates that 
/ ^ and hence J ^ 0. The calculation of I includes only the contribution 
of the fields, but clearly this may, for systems involving obstructions, not 
be sufficient. We can consider adding an external contribution due to the 
potential obstruction and see whether this restores J-conservation. How the 
potential obstructions affect, if at all, the orbital angular momentum needs to 
be considered. The symmetric obstructions can be written as a contribution 
to the potential term V(<j)) in terms of Heaviside functions. The potential 
term and this inhomogeneity expressed in terms of Heaviside functions can 
be written as: 

Vfo) = \ls (1 - <Pt) ± ^ (1 - 01) [Q(y + y ) - Q(y - y )} . (26) 

Using the above definition of V(4>), one can compute the contribution made to 
I in addition to the fields already computed from (J6j). One needs to construct 
q from its constituent parts as shown in (Q: 

' 5W 

Using the properties of the Heaviside functions and their relations to the 
5-function one can show that the potential obstruction's contribution to q is 
given by: 

q = ±</> 3 <9 x 3 r [8(y + y )-$(y- Vo))] ■ 

With this expression for q, we can evaluate the total rate of change of the 
orbital angular momentum due to the obstruction. We find: 

1 = \ J J {x 2 + V 2 ) Q dx dy 

j ~p /'oo roo 

= — dx (x 2 + y 2 )d x (-$){5(y + y )-5(y-y ))}dy 



q = -eijdidiafl = -eijdi — ■ dji 



^ /_~ dx (* 2 + y 2 R(^ 



y=-yo 

(27) 

y=yo 
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The integrals in (127)) have been calculated during each 
simulation. Fig. llTU]) shows the numerically computed integral contributions 
due to the obstruction and the numerically calculated derivative of the orbital 
angular momentum from Fig. (Q plotted with respect to time, for a potential 
barrier with b = 2 and T = 0.25. It can be seen from the plot that the time 
evolution of the integral contributions exactly matches that of I due to the 
fields, so that we have: 



I 



^ fields ^barrier 



d 

~dt 
0. 



1 ±y f°° 1 9 

(x 2 + y 2 ) qdxdy^+ — / dx (x 2 + y 2 ) d x (-(f>l) 



(28) 

y=-yo 

(29) 

y=yo 
(30) 



Due to discretisation effects and numerical inaccuracies, the result is not 
exact but the qualitative features of the integral contributions makes this a 
very consistent result. 




Figure 10: Plots of the numerically calculated time derivative of the total 
angular momentum J = I + m and the contribution of the barrier to / for 
b=2 T = 0.25. 



Fig. (flTl) shows the numerically computed integral contributions due to the 
obstruction and the numerically calculated derivative of the orbital angular 
momentum from Fig. ([8]) plotted with time, for a potential hole with 6 = 2 and 
T = 0.25. It is clear from these plots that the conservation of the total angular 
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momentum J is restored by the introduction of the terms corresponding to 
the potential obstructions' contribution to I. 




Figure 11: Plots of the numerically calculated time derivative of the total 
angular momentum J = I + m for the fields and the contribution of the hole 
to i for b=2 T = -0.25. 
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VIII Conclusions 



Our studies have shown that the scattering of baby skyrmions of our model 
off potential obstructions, for which the dynamics is governed by the Landau- 
Lifshitz equation, exhibits some nontrivial results. 

We have managed to understand quite well the observed scattering prop- 
erties of our skyrmions despite their, at first sight, somewhat non-intuitive 
behaviour. Thus, in the case of a potential hole the skyrmions were unable to 
penetrate it and so moved parallel to the 'x'-axis at a distance y s max from the 
hole. The energy considerations have shown that the skyrmions, in systems 
involving a potential hole, were no longer bound and so could and did move 
away from each other. In the barrier systems the skyrmions were able to 
traverse the barrier. Our simulations have shown that as the skyrmions tra- 
versed the barrier their distance of separation d decreased to overcome this 
increase in potential energy. At the same time the skyrmions sped up as they 
climbed the barrier. Whether the skyrmions were bound or not could only 
be determined by the details of the energetics of the system. It was found 
that the skyrmions were bound for certain low values of the parameters b 
and T. At higher values they were no longer bound and were free to move 
away from each other. 

An interesting observation of our simulations was the apparent non- 
conservation of the total angular momentum J (given its usual definition). 
This non-conservation of J was due to the non- conservation of the orbital 
angular momentum I, as we have found that in all of the simulations the 
total magnetization in the third direction m was well conserved in time. At 
the same time we showed that I 7^ 0. Thinking about this further we showed 
that when a system possesses potential obstructions these obstructions made 
a significant contribution to I. Hence one has to modify the conventional 
definition of I. We have found this missing contribution and we have shown 
that its change compensates /, resulting in the overall conservation of I and J 
for the full system. We believe that most of the results presented here form 
a generic basis for the description of the scattering of baby-skyrmion con- 
figurations in Landau-Lifshitz models. This is primarily due to conservation 
laws of the Landau Lifshitz systems, as constructed by Papanicolaou and 
Tomaras[10], and the observation that the potential obstructions contribute 
to the conservation laws of the system. 



26 



References 



[1] Brand J, Piette B and Zakrzewski W 2005, Scattering of topological 
solitons on holes and barriers, J. Phys. A 38, 10403-10412 

[2] J.M. Speight 2006, Sigma models on curved space and bubble refraction 
in doped antiferromagnets, Nonlinearity 19 1565-1579 

[3] Bobeck A.H and Delia Torre 1975 Magnetic bubbles (Selected topics in 
Solid State Physics, North-Holland publishing company) 

[4] O'Dell T.H 1974 Magnetic bubbles (The Macmillan Press Ltd.) 

[5] Derrick G.H 1964 Comments on nonlinear wave equations as models for 
elementary particles J. Maths. Phys 5 1252 

[6] Weidig T 1999 The baby Skyrme models and their multi-skyrmions Non- 
linearity 12 1489-1503 

[7] Papanicolaou N and Zakrzewski W 1996, Dynamics of magnetic bubbles 
in skyrme model Phys. Lett. A 210 328-336 

[8] Papanicolaou N and Zakrzewski W 1995, Dynamics of interacting mag- 
netic vortices in a model Landau-Lifshitz equation Physica D 80 225-245 

[9] Eslami P, Sarbishaei M and Zakrzewski W 2000, Baby Skyrme models 
for a class of potentials Nonlinearity 13 1867-1881 

[10] Papanicolaou N and Tomaras T.N 1991, Dynamics of magnetic vortices 
Nuc.Phys B 360 425-462 



27 



